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The theory of electron-phonon interaction in the presence of strong correlation has been inves- 
tigated in the present work. Due to the so called spin-charge separation, it is argued that the 
electron-phonon interaction in the strongly interacting systems is fundamentally different from the 
weakly interacting counterparts. Both spinous and holons are shown to carry electric charge and 
couple to the electromagnetic field and consequently to phonons. A systematic approach to calculate 
the dynamical charge of spinous and holons is presented, followed by a discussion of implications of 
this model for the high temperature superconductors. It is demonstrated that the unusual oxygen 
isotope effect in high Tc materials can be easily understood within this framework. 
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I. INTRODUCTION 

Providing a satisfactory description of the strongly cor- 
relating systems is one of the most challenging problems 
in theoretical condensed matter physics. The discovery of 
the high temperature superconductivity in cuprate com- 
pounds^'^ has made this challenge even more intriguing. 
Although, there has been a huge progress in this field 
during the past three decades, we still lack enough knowl- 
edge about the nature of the electron-phonon interaction 
in the strongly correlating systems. In the present paper, 
this issue is addressed from the spin-charge separation 
point of view and the effect of strong interaction on the 
coupling of electrons to phonons is investigated. 

Although at the moment there is no general consensus 
about the underlying mechanism of superconductivity in 
cuprates and the glue that pairs up electrons, there are 
several theories that more or less capture the essence of 
physics in these materials. However, it should be men- 
tioned that it is widely accepted that the physics of high 
Tc is that of doped Mott insulator. In this paper, we 
particularly work in the framework initially introduced 
by P. W. Anderson. According to the Anderson's idea 
of preformed Cooper pairs, the basic physics of high Tc 
cuprates is governed by strong correlation physics due 
to their proximity to the Mott insulator phase at half 
filling^. This approach is based on the resonating va- 
lence bond(RVB) state and the spin charge separation 
picture. At half filling, RVB state is made up of singlet 
Cooper pairs, but due to strong correlations, they do not 
have phase coherence and therefore, we obtain an insu- 
lating phase. Upon doping the system, they acquire the 
required phase coherence and we obtain a superconduct- 
ing state. 

This approach can explain basic properties of cuprates 
and yields a phase diagram similar to that of high Tc 
cuprates. On the other hand, observation of strong iso- 
tope effect on both the transition temperature and the su- 
perfluid density in cuprates'*"^" , indicates the importance 
of the electron-phonon interaction in understanding the 
physics of high Tc^^"^'^. The reported oxygen isotope 




FIG. 1. (Color online) Schematic phase diagram and the elec- 
tric charge of holons and spinous. SM stands for the strange 
metal phase where neither holons nor spinous condense. Sc is 
superconducting, FL fermi liquid and PC pseudogap phase. 



effect (OIE) is very different from that of conventional 
BCS superconductors in which the electron-phonon in- 
teraction is responsible for the pairing mechanism. These 
observations led us to reexamine the electron-phonon in- 
teraction in the presence of strong correlations and its 
effect on the pairing mechanism. 

The deformation of a lattice increases its electric po- 
tential energy quadratically in the simplest approxima- 
tion. Phonons are the quantum of these vibrations. A 
charged quasiparticle e.g. an electron couples to the 
phonon field proportional to its electric charge. There- 
fore, electrically neutral quasiparticles do not couple to 
phonons. A fundamental question about the electron- 
phonon interaction in the strongly correlated materials 
that was raised by the spin charge separation assumption 
is how to assign electric charge to spinous and holons. Do 
spinous carry all the electric charge and holons remain 
neutral? Should we assume the electric charge goes to 
holons and spinous are neutral quasiparticles? The main 
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goal of this paper is to find the answer to this question 
unambiguously. The only assumption made here is the 
spin-charge separation and we demonstrate that the elec- 
tric charge is a dynamical quantity that changes around 
the phase diagram of high Tc compounds. For exam- 
ple, in the Fermi liquid phase, we show that only spinous 
carry electric charge while in the pseudogap phase only 
holons carry the electric charge. In the superconducting 
region both holons and spinon carry an electric charge. 
In every region of the phase diagram qh — Qs = +1) where 
Qh is the electric charge of holons and Qs is that of spinous. 



II. METHOD 

In the slave boson method, electron operator at site i 
with spin a is written in terms of slave particles as 



= h- 



(1) 



Holons are assumed to be hard-core bosons. Accord- 
ing, there is no doublon state in the Hilbert space and at 
each site we either have spinon or holon. In other words 



1. 



(2) 



The above definition for the electron operator is clearly 
invariant under local U(l) gauge transformation, pro- 
vided spinons and holons carry the same charge with 
respect to the internal gauge field. By scaling the in- 
ternal gauge field we can assume the internal charge of 
slave particles is eint = 1 (in units of electron charge e). 
Due to the local constraint on the Hilbert space, hopping 
of a spinon from site i to site j, should be accompanied 
by another hopping of a holons from site j to site i. So 
we conclude that the current carried by spinons is equal 
but opposite to the current carried by holons and they 
add up to zero. Consequently, we have 



-J, (x) = 0. 



(3) 



Now let us assign electric charges qh and qf to holons 
and spinons respectively. On the other hand the electron 

operator transforms under electromagnetic gauge trans- 
formation as ct^ — >■ exp (— iaj) cj^. But we have 



4,a = fla^i ^exp {i {qf -qh)ai) fl^h, 
= exp(i {qf -qh)ai)cl 



Therefore, we have 



Qh-qj = +1- 



(4) 



(5) 



The above constraint does not specify the values for 
qf and qh independently and states that any value for 
qh and qf are accepted as long as they satisfy the above 
equation. In order to satisfy the constraint in Eq. [3] , the 
external electromagnetic field has to generate an internal 



magnetic field proportional to itself as well. It is this 
internal magnetic field that finally resolves the dilemma 
related to the values of qh and qf and uniquely dictates 
their values. In the following we comment more on the 
procedure leading to this conclusion. 

Let us assume a local electromagnetic field 
{x) and internal gauge field {x). Holons 
couple to {qhCLfj, {x) + (/int.ha^ {x)) and spinons feel 
{qfAf_,{x) + qint,{af_,{x)), where qint,h = qint.i = 1- 
Within linear response theory, they generate current as 



3l {x) - n'' {qhA^ {x) 
{x) = {qfA^ {x) - 



(6) 



where we have assumed —i 
equation [3], we have 



JIJ. TV 

•Jf,h' '^f,h 



= ^f,hSi^,v Using 



qhU^+qfllf 
Substituting the above relation in Eq. [6], 



(7) 



(x) = {qh - qf) 



.n^'A^ {x) 



Jji (x) = - {qh 



(8) 



Now using the fact q/, — qf = 1, the above relations 
mean the physical response of the slave particles to the 
electromagnetic field is independent of the initial as- 
sumption for qh and qf. However, we can unambigu- 
ously assign the following physical charges to spinons and 
holons due to their response to the electromagnetic field 
as Jf'"^ = q}%Uf,hA 



qh 



If 



„eff 

qh 



nf 
1 = 



n" 



(9) 



which are independent of the details of charge assigna- 
tion. From now on we drop eff symbol and simply use 
qh and qf. Since H/^^ = —icoaf^h, we can alternatively 
write 



qh 



„eff 



a 



For the physical electric current we have 



(10) 



^ph 



= qh^h + qfJf 



'Jf = {qh-qf)3h = ^^E{n) 

c^h + o-f 



The above relation is known as loffe-Larkin recombina- 
tion formula. Now let us consider different regions of the 
phase diagram of compute the effective charge of holons 
and spinons (see Fig. 1). 
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A. Fermi Liquid Phase 



In the Fermi liquid (FL) phase, holons are condensed 
but spinous arc not. So wc have (h) ^ and Ay = 0, 
where A/ = {fi,ffj,i) is proportional to the pseudo- 
gap energy. In this case, II/j = superfluid density = 
""'m*'"^^ 7^ ^ ( '^'»>c is the condensation fraction and to^ 
effective holon mass, x doping and T temperature), while 
11/ = —iivaf 0, we have 



T'c,h 



{x,T) 



n, 



<ih 



If 



(12) 



Therefore, in the Fermi liquid wc have g/ ~ —1 and 
therefore, they have nonzero overlap with physical elec- 
trons and carry the same charge. On the other hand 

~ and wc can safely assume that holons are electri- 
cally neutral in the FL phase. This can also be directly 
seen from c| ^ ~ {h) fj ^. One important result of this ar- 
gument is that in the absence of the pscudogap («. e. when 
As = 0), condensed holons do not couple to phonons, 
since phonons only couple to electrically charged quasi- 
particles. In other words, phonons create local electro- 
magnetic field and this field induces another (internal) 
gauge field. Holons couple to the sum of these two fields. 
The induced internal field almost cancels out the electric 
field and therefore, holons do not response to the elec- 
tromagnetic field created by phonons. Therefore, when 
in the FL phase where holons condense and pseudogap 
is absent, {i.e. at low enough temperatures and large 
enough doping) only spinous respond to the phononic 
electromagnetic field. This argument is valid close to 
the phase transition between the Superconducting (SC) 
state and the FL state, i.e. in the overdoped region. In 
the next section we discuss that the phenomenology of 
cuprates inspires us to assume the isotope effect is solely 
due to the holon-phonon interaction. This immediately 
explain why there is no isotope effect on Tc in the over- 
doped region where SC state undergoes transition into 
the FL state. 



B. Superconducting Phase 

In the superconducting (SC) phase , both holons and 
spinous are condensed. So we have = (hh) ^ and 
Ag ^ 0. In this case, since 11 = , both and qf are 
nonzero. On the hand, for the d-wave superconductors 
we have 




(14) 



Computing the effective charge of spinous and holons 
and we obtain qf oc —lih oc —x and qh, ocUf oc {1 — x) 
respectively. Therefore, in the superconducting state, 
spinons only respond to the x fraction of the electromag- 
netic field. Since holons are effectively charged quasipar- 
ticles in this state, (or since the internal gauge field does 
not cancel out the local electromagnetic field of phonons 
completely) they couple to phonons. The larger pseudo- 
gap value, the stronger interaction between holons and 
phonons. This is an intuitive way to justify why the 
isotope effect is directly dependent to some power of the 
pseudogap energy in the experimental observations. This 
means that gauge fluctuations, renormalizes the coupling 
constant of holon-phonon interaction from the brae value 
7, to 7* = A"7. Using this expression and since A^ de- 
creases with doping, we can easily explain why the iso- 
tope effect is a decreasing function of doping. Moreover, 
in the overdoped region, Ag = at the transition temper- 
ature and therefore, holons do not interact with phonons. 
Subsequently, we do not expect isotope effect on the T^. 
in the interface between SC and FL phases. However at 
T = 0, As even in the overdoped region and the mass 
of holons enhances in an isotope dependent way and this 
explains the nonvanishing isotope effect on the superfluid 
density and the London penetration depth in this region. 



C. Pseudogap phase 

In the pseudogap phase (PG), spinons form Cooper 
pairs with d-wave symmetry and condense. This region 
is above the underdoped SC region. In the underdoped 
side of the SC phase, according to phase fluctuation sce- 
nario for the SC transition, superconductivity disappears 
above Tbkt (x) , the temperature above which Cooper 
pairs lose their long range phase coherence and their 
phase fluctuates strongly, vortices proliferate and super- 
conductivity disappears. We assume this temperature is 
of the order of BEC temperature for the holon gas, but 
a bit less i.e. Tbkt ^ 2bec- Therefore, we assume in the 
PG region holons are still condensed between Tbec {x) 
and Tbkt {x). Above Tbkt we have {h) = 0. Accord- 
ingly 

Tbkt<T<Tbec: qh{x,T) = „ < 1, 9/ < 

11/ -I- fl/i 



n,j{x,T) = ncj{x,0)(l-,^'j (x A/ (a;, T)" (13) qUx,T) <x n,j {T,x) <x (^1 - <xA}{T,x) (15) 



with some exponent a 2. At zero temperature, al- 
most all spinons condense and we obtain ricj ~ 1 — x. 
For holons we have 



However, above the BEC temperature, since 
the superfluid density of holon gas vanishes i.e. 

lim^j_).o n (g = 0, w) — >■ 0, we have 



4 



T > T, 



BEC 



1. 



<ih 



0(16) 



Therefore, deep in the PG region, only holons carry 
electric charge and couple to phonons. Spinons carry 
electric charge only in a narrow region between Tbkt {x) 
and Tbec {x). 



III. OXYGEN ISOTOPE EFFECT IN 
CUPRATES 

It has been experimentally verified that ^^O/^^O iso- 
tope substitution affects the SC transition temperature 
as well as the superfluid density. Since oxygen is lighter 
than copper in cuprate compounds and displaces easier, 
oxygen atoms contribute more to the lattice vibrations 
and cause isotope effect. Experimental observations indi- 
cate a strong OIE on Tc only in undcrdoped cuprates(see 
Fig. 1), while there is no considerable OIE in overdoped 
cuprates. On the other hand OIE on the superfluid den- 
sity has been reported for both sides^^""'^^. In both cases 
the isotope exponent decreases as we approach the op- 
timal doping from the underdoped side. On the other 
hand, there is no isotope effect on the effective mass of 
the nodal quasiparticles^^. It is worth comparing these 
behaviors with that of conventional BCS superconduc- 
tors. In conventional superconductors, there is no iso- 
tope effect on the superfluid density as well as the effec- 
tive mass of quasiparticles, and the isotope effect on Tc 
is independent of doping with isotope exponent around 
1/2^^. We would like to mention that the BCS theory 
and its generalizations e.g. Migdal-Eliashberg theory are 
based on the adiabatic approximation. Adiabatic limit 
assumes the phonon energy is much less than the typical 
energy of electrons i.e. E-p. Therefore, adiabatic approx- 
imation in the Einstein model means <C -Ep, where 
LOj^ is energy of Einstein phonon mode (optical phonons) . 

In high Tc cuprates, the situation is different if we 
assume spin charge separation. The typical energy of 
spinons is of the order of electrons. Experimental data 
implies that the kinetic energy of electron (and as a re- 
sult that of spinons) is around the exchange energy J 
which is about 1500 Kelvin. On the other hand, the typ- 
ical energy of phonons is expected to be several hundred 
Kelvin. Therefore, E-p > lu^^ and we are in the regime of 
adiabatic spinon-phonon interaction. Subsequently, the 
mass of spinons renormalizes as ~ (1 -f- A), where 
A is the coupling constant of spinon-phonon interaction. 
In the simplest model it can be shown that A ex , and 

E 

^ A is isotope independent. Therefore, 



smce Wjj, oc 

oxygen isotope substitution does not enhance the effec- 
tive mass of spinons. This is consistent with the exper- 
iment where there is no isotope effect on Fermi velocity 
by Laser ARPES, while there is shift in kink energy 

For holons we must be more careful. Holons are 
bosonic objects and their typical energy is very small 



compared to their bandwidth as long as they condense. 
Therefore, for holons we have Eh << and we should 
use the results of the non-adiabatic limit of the holon- 
phonon interaction. 

To explain the isotope effect in cuprates, note that in 
the slave boson approach, superfluid density is inversely 
proportional to the mass of holons, pa,c{x,T = 0) = 

Since the isotope substitution changes the superfluid den- 
sity in the whole region and decreases with increasing 
doping (it finally vanishes close to the FL phase) , we con- 
clude that the mass of holons should change with isotope 
substitution as well. On the other hand, Tc is propor- 
tional to the superfiuid density in the underdoped region 
and therefore, Tc{x) cx 1/toJ^ in the underdoped side. 
This immediately means we can explain OIE on Tc pro- 
vided we show OIE on m^. On the other hand, Tc in the 
overdoped region is given by the transition temperature 
of spinons above which pseudogap closes. Since antiferro- 
magnetic exchange dominates over the phonon mediated 
attraction in cuprates, oxygen isotope substitution does 
not affect Tc much in the overdoped region. 

From the above discussion we conclude that an iso- 
tope dependent holon mass can easily explain the ob- 
served OIE in cuprates. In the first half of this paper, 
we discussed that phonon field couples to charged quasi- 
particles. Therefore, holons can in principle couple to 
phonon if they acquire nonzero charge. In the previous 
section, we demonstrated that qh{x,T) cx pf{x,T) cc 

^1 — jjl^^ ■ Therefore, we can write down the following 

effective model for holons 



<4J> 



70 



--hl^hn (ajj + a„) 



(17) 



where 70 is the coupling constant of electron-phonon 
interaction and al (a^) is the creation (annihilation) 
operator of phonons at site i. Note that qh fac- 
tor renormalizes bare 70 to qhjo- The above model 
has been extensively studied. In the non-adiabatic 
regime, we can use the results of single polaron 
theory^"'^^. Therefore, we can do the powerful Lang- 
Firsov transformation^^. 



I.e. 



H 



H 



'He'- 



J2q,n IhlQhlK {ai - a„) . Un- 



where: S = 7== , , 

der this transformation, hi transforms as: e~^hie^ 

7o(a|-ai) 

tor by its average, we finally have: 



hiBxp qh 



If we replace the exponential fac- 



-ff ^ (2i^ (cos ka; + cos ky) + pi) h\hk 

k 

+ ^c., (dtd, + l/2) + ... (18) 
where ellipses stand for interaction terms and 
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(19) 
(20) 

(21) 



in which is the oxygen mass. It is clear that Aq is 
mass independent from its definition. Now if we expand 
the energy of holons around fc = we have 



eh (k) 



2ml 



"^h {x, T) = nih exp [g^ (x, T)) , 



(22) 
(23) 



where = As we see from the above equation, 



scattering off the phonons enhances the mass of holons 
by an exponential factor. This factors depends on the 
isotope mass and as a result we expect the OIE on any 
physical quantity depending on the mass of holons. The 



OIE on is 



dlogM^ 



(24) 



/? reflects the same behavior which is consistent with em- 
pirical data. 

B. Oxygen isotope effect on Tc 

The oxygen isotope effect on Tc is determined by the 
a isotope exponent which is defined as: 



a {x) = — 



dlogTcjx) 
d log Mo 



(27) 



In the underdoped region, phase fluctuation studies 
show that Tc is controlled by the superfluid density i.e. 
Tc ^ 2bkt oc Ps cx \/m*j^. Therefore, we have: 



X < xoD : 

^^^^^dlogm*(.,T, 



d log Mo 



Qh (a;, Tc^ 



Ao 



(28) 



where xqd stands for the optimal doping. Since 
qh{x,Tc) (X pf{Tc) is a decreasing function of doping 
and at the optimal doping and beyond vanishes, isotope 
exponent dies off as we approach the optimal doping and 
finally becomes negligible. On the overdoped cuprates, 
at the boundary between the SC and FL phases qh van- 
ishes and therefore, we do not expect isotope effect on Tc 
for that case. These conclusion are also consistent with 
experimental observations. 



A. Oxygen isotope effect on the superfluid density 



IV. SUMMARY AND CONCLUSION 



The oxygen isotope effect on the superfiuid density is 
defined as: 



p{x) 



dlogps {x,T^O) 
dlogM„ 



(25) 



According to the loffe-Larkin formula the physical 
superfluid density is related to the superfluid density of 
spinous and holons as: p~i^ — pf^^ + pj^ . Since conden- 
sation fraction of holons and spinous at zero temperature 
are x and 1 — x respectively, and from p — ^ we have: 

Pph i^) ~ ^ + I^- ^'^^ small values of x, we have 
Pph (0) ~ Therefore, we have: 



Pix) 



dlogm^ (a;,0) 
d log Mo 



Ao 



(26) 



Since qh(x,T) cc llf{x,0) — Uc.f {x,Q) /m*^ and the 
condensation fraction of spinons decreases with increas- 
ing doping until vanishes as we approach the FL phase. 



We have studied the electron-phonon interaction in the 
strong correlation regime within the spin charge separa- 
tion framework. It was shown that by the careful study 
of the gauge theory and by taking local constraints into 
account, the electric charge of slave particles can be de- 
termined unambiguously. The charge of holons was cal- 
culated and shown to be proportional to the superfluid 
density of spinons in the SC region. Furthermore, the 
charge vanishes in the FL region whereas it grows in the 
PG phase. We argued that holons couple to phonons 
due to their nonzero electric charge. The holon-phonon 
interaction should be treated in the non-adiabatic limit 
since holons are bosonic quasiparticles and their energy is 
much less than that of phonons as long as they condense. 
The implications of our arguments on the OIE in high Tc 
cuprates were discussed above, and we have shown that 
the effective charge of holons and the exponential depen- 
dence of their mass on qh readily explains the unusual 
16q^18 isotope effect in cuprates. Our model successfully 
explains the general trend in experiments and supports 
the spin charge separation idea. 
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